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ABSTRACT 

Context. The winter seeing at Concordia is essentially bimodal, excellent or quite poor, with relative proportions that depend on 
altitude above the snow surface. This paper studies the temporal behavior of the good seeing sequences. 

Aims. An efficient exploitation of extremely good seeing with an adaptive optics system needs long integrations. It is then important 
to explore the temporal distribution of the fraction of time providing excellent seeing. 

Methods. Temporal windows of good seeing are created by a simple binary process. Good or bad. Their autocorrelations are corrected 
for those of the existing data sets, since these are not continuous, being often interrupted by technical problems in addition to the 
adverse weather gaps. At the end these corrected autocorrelations provide the typical duration of good seeing sequences. This study 
has to be a little detailed as its results depend on the season, summer or winter 

Results. Using a threshold of 0.5 arcsec to define the "good seeing", three characteristic numbers are found to describe the temporal 
evolution of the good seeing windows. The first number is the mean duration of an uninterrupted good seeing sequence: it is tq = 
7.5 hours at 8 m above the ground (15 hours at 20 m). These sequences are randomly distributed in time, with a negative exponential 
law of damping time ti = 29 hours (at elevation 8 m and 20 m). The third number is the mean time between two 29 hours episodes. 
It is r = 10 days at 8 m high (5 days at 20 m). 

Conclusions. There is certainly no other site on Earth, except for the few other high altitude Domes on the Antarctic plateau, that can 
get close to these really peculiar seeing conditions. 
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1. Introduction 

Regarded as astronomical sites, the highest points of the 
Antarctica plateau present many obvious advantages due to the 
local climate and the remoteness from any polluting civilization. 
They also benefit from an interestingly unique distribution of tur- 
bulence. This has been extensively measured at Dome C since 
the first winter-over permitted in 2005 by the French - Italian 
Concordia station operation (Aristidi et al, 2009). Winter and 
Summer display very different but both unusual vertical distribu- 
tions of the turbulent energy. Generally speaking, the situation is 
dominated by the presence of a surface inversion layer that be- 
comes very turbulent when the temperature gradient is strong in 
winter, and can completely vanish in summer when this gradi- 
ent becomes flat. In summer it depends on the Sun's elevation, 
and is then strongly time dependent, with an optimum period of 
a few hours of excellent seeing every day in the middle of local 
afternoon (Aristidi et al, 2005). In the other 3 seasons, the mean 
seeing is essentially altitude dependent above the snow surface. 
Above this surface turbulent layer that contains, statistically, 95 
percent of the total Cn^ along the line of sight, it has proved to be 
statistically independent of the season, within the measurements 
accuracy. Its mean value is between 0.3 and 0.4 arcsec as soon as 
the telescope is located above a sharply defined altitude thresh- 
old, that fluctuates around a mean value of the order of 25 m. The 
consequence is that the non summer seeing displays a nearly bi- 
modal statistical distribution. It is indeed as good as 0.3 to 0.4 
arcsec 50 percent of the time at 25m above the surface, this frac- 
tion of time decreasing to about 40 percent at 20m and slighdy 
less than 20 percent at 8m. But it is obviously not equivalent to 
have 40 percent of good seeing spread in many short sequences 



of seconds to minutes rather than distributed in extended long 
sequences measured in hours or days. This paper addresses the 
temporal distribution of this good seeing percentage. It goes be- 
yond the first analysis made in Aristidi et al. (2009), using a 
method to compensate for the gaps in the data. 



2. The autocorrelation method 

This paper exploits the Concordia DIMM data. A DIMM located 
at about 8m above the snow surface is permanently operated 
since the end of 2004, after having already provided a summer 
season one year before. A second DIMM was set at 20m high on 
the roof of the Concordia station during 3 months in 2005. There 
is also a GSM, made of two DIMMs on the snow surface. They 
are not exploited in the present paper. 

The DIMM seeing data sets have a sampling rate of 2 min- 
utes. Their data time series are not continuous, being relatively 
often interrupted either by adverse weather, or by various tech- 
nical reasons such as frost on the optics, electronics or computer 
shut downs, loss of star tracking, .... It is then difficult to directly 
track the continuity of good seeing sequences in the original data 
files as a gap in such a sequence of good seeing data can be due 
either to a bad seeing moment or to a simple lack of data. This 
study must then be made on a statistical basis. 

The autocorrelation method applied to data windows proved 
to be a very efficient tool for that purpose. Two different temporal 
windows can be defined, by means of a function equal to 1 or 0. 

The first one is called vvj.(t) and defines the existence of data 
at a time t. Weit)-l if data exist, otherwise. 
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The second one Wg(t) defines the seeing quaHty. It uses a 
threshold cTg so that Wg(t) =1 if the seeing cr(f) < o-g and Wi;(t)=0 
if cr(f) > cr^. In all the rest of this paper, we will use o-g = 0.5 arc- 
sec, a value that is near the minimum of the gap between good 
and bad seeing in the winter seeing histograms (Aristidi et al, 
2009). 

This second window function, however, is not directly ac- 
cessible since a significant fraction of measurements is missing. 
The only accessible one is the product of these two windows, 
that will be 1 when the measurement does exist and provides a 
good seeing value, and when either the seeing is not good or 
the measurement does not exist. This accessible window func- 
tion can be called Wge{t)-We{t) Wg{t) 

It is unfortunately not possible to recover the interesting but 
unknown window Wg{t) as a division by We{t) that contains many 
zeroes is not feasible. Such a division, however, is possible in the 
domain of the autocorrelations, as long as the data sets are long 
enough regarding their characteristic time. Assuming that Wg and 
We are statistically independent, the autocorrelation of Wge{t) can 
be written as 

rgfW ^ j;^^E[wMWg(t)We(t + T)Wg(t + T)] (1) 

= ^^EM) Weit + t)] E{wg{t) Wg(t + t)] 

= r,(T)r,(T) 

where TgiT), Tgir) and r^p(T) are the autocorrelation functions 
of We(t), Wg(t) and Wgdi) respectively, irig and nig are the mean 
values of We and Wg (the mean value is here equal to the vari- 
ance since the functions and Wg equal or 1) and E[] denotes 
the expected value. This technique has been extensively used in 
the past to analyse heliosismic data (Lazrek, 1993) and is in fact 
a method of deconvolution in the Fourier plane. As we benefit 
of very long time series, spread over several years of data that 
can be cumulated, and as there is very little suspicion of depen- 
dence between the observation window w^it) and the good see- 
ing occurrences Wg(t), we can use this division and thus study 
the statistical properties of the good seeing sequences through 
the autocorrelation Fgir) of the corresponding window Wg{t). 

It is now well established that the summer and winter condi- 
tions at Concordia are drastically different. We are mostly inter- 
ested in this paper by the winter conditions. The "winter" in this 
study will be defined as the 6 months during which no significant 
temperature variation is visible, which is from early April to the 
end of September We will also qualify and illustrate the method 
by a summer analysis, the summer data being taken during the 
permanent sunlight season, from early November to February 
10th. 

Figures 1 and 2 illustrate the data window function autocor- 
relation TeiT). During the first two winters, in 2005 and 2006, 
only one astronomer was present on the site. It is clear on Fig. 
1 that the automation of the instrument was not fully successful 
yet, as demonstrated by the one-day periodicity of the window 
function. In 2007 for the first time, two astronomers wintered 
over and could then take care of the instrument more perma- 
nently. The one-day periodicity of the window function is not 
totally cancelled, but almost. 

The overall filling factors of the data sets are provided by the 
asymptotic values of the autocorrelations functions. It is indeed 
well known (Papoulis, 1984) that r^iT) — > Elw^it)^ Im^ - 
where t — > oo, as the quantities We{t) and We{t + t) become 
statistically independent. The autocorrelation tends towards the 
mean value of Wg, i.e. the time percentage where Wg - \. This 
filling factor was found to be about 40 percent during the first 
two years and slightly more than 30 percent in 2007, due to a 
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Fig. 1. Autocorrelation F^ (r) of the existing data window, aver- 
aged on the 2005 and 2006 winter seasons. The one-day period- 
icity is clearly visible. 
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Fig. 2. Autocorrelation T^iT) of the winter 2007 data window 
function. A long gap due to an important technical problem oc- 
curred during that winter, so that the overall filling factor is a 
little less, but with two astronomers wintering over instead of 
only one, the amplitude of the 24-hour periodicity is drastically 
reduced. 

relatively long technical interruption during the first half of that 
winter season. 

3. The summer situation 

In summer, there are generally a few more observers, but there 
is also more technical maintenance, so that the gaps in the data 
are not less than in winter The overall filling factors are still in 
the range 30 to 40 percent but with several people taking care, 
the one-day periodicity is almost absent in the window We(t) (not 
shown here). However, the seeing itself has a 24-hour periodicity 
as it displays a minimum every day in the afternoon (Aristidi et 
al, 2005), so that the autocorrelation of the good seeing window 
function Tg{T), obtained by the division of Vgg{T) by Vgij), is 
then expected to display a large amplitude one-day periodicity. 
That is confirmed by Fig. 3, averaged on all available summer 
seasons from 2003/2004 up to 2007/2008 and then statistically 
very robust. 

This figure deserves some comments. The first comment is 
about the asymptotically behaviour of the autocorrelation, which 
oscillates between 0.42 and 0.72, around a mean value of 0.57. 
It means that 57 percent of the time during the 3-month summer 
season the seeing is better than our threshold cTg - 0.5 arcsec. 
Another comment concerns the first 6 hours, where the auto- 
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Fig. 3. Autocorrelation rg(T) of the good seeing data window, 
averaged on all the summer seasons from 2003/2004 up to 
2007/2008. This time, the one-day periodicity is not due data 
gaps but to a real daily variability of the seeing. 
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Fig. 4. Autocorrelation Vgir) of the good seeing windows, av- 
eraged on 4 winter seasons seasons from 2005 to 2008. The 5 
parameters fit described by Eq. 2 is superimposed to the curve. 



correlation shows an almost exactly linear decrease. Such a lin- 
ear decrease corresponds to a rectangular window function. It 
thus means that a significant part of the summer good seeing se- 
quences is indeed made of continuous sequences of 6 hours long 
on average, every day at the same time. This is exactly what we 
already know of the summer seeing (see Fig. 3 in Aristidi et al, 
2005), which has been measured to be better than 0.5 arcsec be- 
tween 2 and 8 p.m. every day in summer 

It can be noted that 6 hours per day represents 25 percent, 
while the total fraction of good seeing data is statistically more 
than twice this percentage. There are then many other episodes 
of good seeing in summer They are shorter and more randomly 
distributed. For instance, the quick decrease of the autocorrela- 
tion from 1 to ~ 0.9 in the first few minutes indicates that an 
order of 10 percent of good seeing or bad seeing events occur as 
isolated events of one or a very few consecutive 2-minute mea- 
surements. 

Besides that last remark, this autocorrelation analysis of the 
summer data does not provide many new informations indeed. 
However it does qualify the method and will be used as a refer- 
ence for the winter data understanding. 



4. The winter situation 

Let's start this section by looking at Fig. 4, the mean autocorre- 
lation function r^(T) of the good seeing windows, averaged on 4 
winter seasons from 2005 to 2008. This is the seeing measured 
on the Concordiastro platform at 8 m. Thanks to the 24 months 
of data exploited here, this autocorrelation is statistically very 
robust again. 

The asymptotic value, slightly less than 0.2, is a good ap- 
proximation of the probability of obtaining good seeing values, 
i.e. being above the turbulent boundary layer at 8 m high. It con- 
firms the estimation of 18% made by means of the histogram 
integrals (Aristidi et al, 2009). 

The behavior at the origin shows a quick drop from 1 to about 
0.8. It confirms the existence of very short sequences (lasting 
typically less than 10 mn) where the seeing is either continu- 
ously good or bad. These individual events were noticed by the 
observers, and represent 20% of the data. Apart these first points, 
the graph of Tgir) shows three main features. A linear decrease 
in the first 7 hours, a negative exponential distribution of charac- 
teristic time of a few tens of hours, and an horizontal asymptot 



for large values of t. We adjusted to T'gir) the following 5 pa- 
rameters fit : 

^^■«-^(i-^)n(^^)-^^"^-^ (2) 

The first term accounts for the linear part and is truncated by 
the rectangle function H (57;;) equalling to when t > tq. We 
found: 

- The asympotic value A - 0.16. It represents the overall prob- 
ability to observe a good seeing 

- B - 0.52 and ti = 29 hours, the parameters of the exponen- 
tial decrease 

- C = 0.08, To = 7.5 hours, the parameters of the linear de- 
crease 

The linear part is similar to what is obtained in summer It is 
indeed the autocorrelation of a rectangular window function. It 
shows that many good seeing sequences occur in continuous 
runs lasting tq = 7.5 hours on average. An important difference, 
though, is that in summer, such sequences of 6 hours reproduce 
every day at the same time, while in winter, there no signifi- 
cant daily periodicity. The exponential decrease that follows this 
initial linear behavior is the typical winter feature and must be 
interpreted. 

We propose the following model for the temporal window 
function Wg{t). This model described hereafter depends on sev- 
eral parameters which were adujsted so that its autocorrelation 
reproduces correctly the autocorrelation Cg(T) of the data. Wg(T) 
is modeled by well separated episodes including several rectan- 
gular functions of width 7.5 hours. The mean delay between two 
consecutive episodes is noted T, with T = 10 days. 

Inside an episode the rectangular functions start at times f,, 
where f, is a random variable with a negative exponential distri- 
bution of damping time ri = 29 hours. That means that most 
of the windows appear at the beginning of the episode and tend 
to disappear after a delay of 29 hours from the beginning of the 
episode. The mean number A^^ of rectangular windows inside a 
10 days episode is adjusted to account for the value of B in the 
function Cg(T) (eq. 2). We find Nr - 11. These 11 windows are 
mostly concentrated in the first 29 hours of the episode and there 
can be large overlaps. In our simulation, more than 75 % of the 
good seeing values are concentrated in the first 48 hours of a 10- 
day episode. Fig. 5 shows an example of such a 2-day episode. 
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Fig. 5. A nice sequence of good seeing data obtained in winter 
2005 at 8 m high. During 50 hours, there are only a few nearly 
individual values above the 0.5 arc-sec threshold, but also an 11- 
hour gap. 



In addition to these episodes, and to account for the initial 
quick drop of TjCt), we simply simulated the initial drop by 
adding a random number A^, of isolated either good or bad seeing 
(lasting two minutes). A^,- is adjusted to 0.22. 

The general figure at 8 m is then that every 10 days on av- 
erage, there is a two to three day episode during which many 
nearly uninterrupted runs of good seeing occur with individual 
durations of 7.5 hours. 

The seeing a 20 m high. Only a little more than 3 months of 
data have been collected at this height so that the statistical ro- 
bustness is not as good as it is at 8 m. Fig. 6 shows the autocor- 
relation function rg(r) of the good seeing windows at 20 m. This 
curve is clearly not as smooth as the other one, due to its weaker 
statistical robustness. The one-day periodicity is visible because 
a large part of this data set comes from Spring (from July, 23rd 
to October, 31st) when the Sun is present, only part time but 
a longer and longer part till the end of the run. However, the 
same kind of fit can be attempted and works well. Interestingly 
enough, the 29-hour damping time of the exponential decrease 
is not modified by this higher altitude. That may presumably be 
explained by the fact that it depends mostly on the meteorologi- 
cal general situation, that should not be affected by a few meters 
difference of altitude. On the other hand, the nearly linear de- 
crease that was explained by uninterrupted runs of 7.5 hours is 
still present, with now a twice longer duration of 15 hours. The 
numbers A and B are fitted to 0.36 and 0.27 respectively. The 
value A = 0.36 is coherent with the good seeing fraction of 40 % 
derived from the histograms (Aristidi et al., 2009). In our simu- 
lations, A^,. must be adjusted between 6 and 7. This again implies 
lots of overlaps and a mean fraction of good seeing of about 
90 percent during the first 48 hours of a simulated episode. The 
delay T between consecutive episodes is found twice shorter, 
around 5 days. 

The two delays between consecutive episodes, 10 and 5 days, 
imply that during the 100-day winter time (without any sunrise), 
the mean numbers of such good seeing episodes must be about 
10 at 8 m and 20 at 20 m. 

Statistically, This difference can be explained by assuming 
that half the good seeing episodes at 20 m correspond to a bound- 
ary layer upper limit that goes down below 20 m but not as low 
as 8 m, the other half including situations when it will also spend 
some time below 8 m. The same 29-hour exponential decrease 
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Fig. 6. Autocorrelation Tgir) of the good seeing data window, 
at 20-m high during 3 months in Winter/Spring 2005, with the 
same kind of fit. 

at both altitudes gives credit to our model and its meteorological 
interpretation. 

It would not be realistic to attempt more interpretation of 
that curve, that has to be taken as providing an estimate of what 
may happen at 20 m high. However, that estimate is important 
as one may consider that a possible optimal choice for a future 
high angular resolution instrument at Dome C could be to set 
it at ~ 20 m, together with a well qualified GLAO system. It 
would then benefit of a statistics of free atmosphere seeing of 
about 40 percent of the time during all the dark and cold period, 
concentrated inside episodes of at least two days with several un- 
interrupted runs of individual durations being at least 15 hours. 
A number of the order of 20 such episodes can be estimated to 
occur during all the dark period. That will definitely permit long 
integrations in very good seeing conditions. 

5. Discussion and conclusions 

Thanks to the very long data sets exploited here, the validity 
of the autocorrelation division for correcting the data interrup- 
tions and the final statistical robustness are well validated by the 
summer season analysis and results. At least at 8-m high. The 
reduced amount of data available at 20-m high makes the statis- 
tical robustness of the numbers visibly weaker, but the general 
tendency tends to reinforce the interpretation of the 8-m data. 
The 29-hour exponential decay of the distribution of the good 
runs starting times, assumed to depend on the meteorological 
situation, is found to be the same at both altitudes. The second 
characteristic time tq, that can be regarded as the minimum du- 
ration of nearly uninterrupted good runs, is found twice longer 
at 20 m, i.e. 15 hours versus 7.5 hours at 8 m. Finally, the num- 
ber of episodes of excellent seeing is estimated to be twice more 
frequent at 20m, 20 times versus 10 times per winter, and that 
is also easy to understand with a reasonable assumption on the 
vertical motions of the boundary layer upper limit. 
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